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ABSTRACT 

The "P32" Astronomical Observation Template (ACT) 
provided a means to map large areas of sky (up to 45x45 
arcmin) in the FIR at high redundancy and with sampling 
close to the Nyquist limit using the ISOPHOT ClOO (3x3) 
and C200 (2x2) detector arrays on board ISO. However, 
the transient response behaviour of the Ga:Ge detectors, 
if uncorrected, can lead to severe systematic photomet- 
ric errors and distortions of source morphology on maps. 
Here we describe the basic concepts of an algorithm which 
can successfully correct for transient response artifacts in 
P32 observations. Examples are given demonstrating the 
photometric and imaging performance of ISOPHOT P32 
observations corrected using the algorithm. 

Key words: ISO 



1. INTRODUCTION 

From the point of view of signal processing and photom- 
etry diffraction-limited mapping in the FIR with cryo- 
genic space observatories equipped with photoconductor 
detectors poses a particular challenge. In this wavelength 
regime the number of pixels in detector arrays is limited 
in comparison with that in mid- and near-IR detectors. 
This means that more repointings needed to map struc- 
tures spanning a given number of resolution elements. Due 
to the logistical constraints imposed by the limited opera- 
tional lifetime of a cryogenic mission, this inevitably leads 
to the problem that the time scale for modulation of il- 
lumination on the detector pixels becomes smaller than 
the characteristic transient response timescale of the de- 
tectors to steps in illumination. The latter timescale can 
reach minutes. 

Unless corrected for, the transient response behaviour 
of the detectors will lead to distortions in images, as well as 
to systematic errors in the photometry of discrete sources 
appearing on the maps. In general, these artifacts become 
more severe and more difficult to correct for at fainter lev- 
els of illumination, since the transient response timescales 
increase with decreasing illumination. Compared to the 
IRAS detectors, the ISOPHOT-C detectors (Lemke et al. 
1996) on board the Infrared Space Observatory (ISO; Kessler 
et al. 1996) had relatively small pixels designed to provide 



near diffraction limiting imaging. ISOPHOT thus gener- 
ally encountered larger contrasts in illumination between 
source and background than IRAS did, making the arti- 
facts from the transient response more prominent, partic- 
ularly for fields with faint backgrounds. 

A further difficulty specific to mapping in the FIR with 
ISO was that, unlike IRAS, the satellite had no possi- 
bility to cover a target field in a controlled raster slew 
mode. This limited the field size that could be mapped 
using the spacecraft raster pointing mode alone, since the 
minimum time interval between the satellite fine point- 
ings used in this mode was around 8 s. This often greatly 
exceeded the nominal exposure time needed to reach a 
required level of sensitivity (or even for many fields the 
confusion limit). Furthermore, the angular sampling and 
redundancy achievable using the fine pointing mode in 
the available time was often quite limited, so that com- 
promises sometimes had to be made to adequately extend 
the map onto the background. 

A specific operational mode for ISO - the "P32" Astro- 
nomical Observation Template (AOT) - was developed for 
the ISOPHOT instrument to alleviate these effects (Hein- 
richsen et al. 1997). This mode employed a combination 
of standard spacecraft repointings and rapid oversampled 
scans using the focal plane chopper. The technique could 
achieve a Nyquist sampling on map areas of sky ranging up 
to 45x45 arcmin in extent (ca. 70x70 FWHM resolution 
elements) on timescales of no more than a few hours. In 
addition to mapping large sources, the P32 AOT was ex- 
tensively used to observe very faint compact sources where 
the improved sky sampling and redundancy alleviated the 
effects of confusion and glitching. 

In all, over 6% of the observing time of ISO was de- 
voted to P32 observations during the 1995-1998 mission, 
but the mode could not until now be fully exploited sci- 
entifically due to the lack of a means of correcting for the 
complex non-linear response behaviour of the Ge:Ga de- 
tectors. Here we describe the basic concept of a new algo- 
rithm which can successfully correct for the transient re- 
sponse artifacts in P32 observations. This algorithm forms 
the kernel of the "P32TOOLS" package, which is now pub- 
lically available as part of the ISOPHOT Interactive Anal- 
ysis package PIA (Gabriel et. al 1997; Gabriel & Acosta- 
Pulido 1999). Information on the algorithm, as well as the 
first scientific applications, can also be found in Tuffs et 
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al. (2002). The user interface of P32TOOLS is described 
by Lu et al. (this volume). After a brief overview of rel- 
evant aspects of the P32 AOT in Sect. 2, we describe 
the semi-empirical model used to reproduce the transient 
response behaviour of the PHT-C detectors in Sect. 3. 
Sect. 4 describes the algorithms used to correct data. Ex- 
amples demonstrating the photometric and imaging per- 
formance of ISOPHOT P32 observations are given in Sect. 5, 
based on maps corrected using P32TOOLS. 

2. THE P32 AOT 

The basic concept of the P32 AOT, in which imaging was 
achieved using a combination of standard spacecraft re- 
pointings and rapid oversampled scans using ISOPHOT's 
focal plane chopper, is summarised in the contribution by 
C. Gabriel in this volume. There were two basic require- 
ments, outlined in the original proposal for the AOT (Tuffs 
& Chini 1990) which led to the final design: 

The first goal was to give ISOPHOT the capability of 
achieving a Nyquist sky sampling (AO = A/2D = ITarcsec 
at 100 /xm) on extended sources subtending up to 50 res- 
olution elements in a tractable observation time. As well 
as providing a unique representation of any arbitrary sky 
brightnesss distribution, this also improved the effective 
confusion limit in deep observations, which, in directions 
of low background was ^ 100 and 20mJy rms at 200 and 
100 Atm, respectively, for ISOPHOT-C. 

The second requirement was to achieve a redundancy 
in the observed data. This led to the sampling of each 
of a given set of sky directions by each detector pixel at 
different times. The redundancy was achieved in part by 
making several (at least 4) chopper sweeps at each space- 
craft pointing direction. Further redundancy was achieved 
by having an overlap in sky coverage of at least half the 
chopper sweep amplitude between chopper sweeps made 
at successive spacecraft pointing directions in the Y space- 
craft coordinate. Thus, a given sky direction was sampled 
on three different time scales by each detector pixel: on 
intervals of the detector non-destructive read interval, on 
intervals of the chopper sweep, and on intervals of the 
spacecraft raster pointings. Although the redundancy re- 
quirement was originally made to help follow the long 
term trends in detector responsivity during an observa- 
tion, in practice it also proved very useful in removing ef- 
fects caused by cosmic particle impacts (Gabriel & Acosta- 
Pulido 2000), the so-called "glitches" , from the data. This 
was especially important to optimise the sensitivity of 
deep observations using the CI 00 detector, which are gen- 
erally limited by glitching rather than by confusion. 

Both these requirements led to dwell times per chopper 
plateau which could be as low as 0.1s. For a given sky 
brightness, this also had the consequence that the detector 
was generally read faster than in the other AOTs using the 
PHT-C detectors. 



2.1. The P32 "natural grid" 

The fact that the spacecraft pointing increment in Y was 
constrained to be a multiple of the chopper step interval 
means that a rectangular grid of directions on the sky can 
be defined such that all data taken during P32 measure- 
ment fine pointings will, to within the precision of the 
fine pointings, exactly fall onto the directions of the grid. 
This so-called P32 "natural grid" allows maps to be made 
without any gridding function, thus optimising angular 
resolution. The P32 natural grid is also a basic building 
block of the algorithm to correct for the transient response 
behaviour of the detector, as the algorithm solves for the 
intrinsic sky brightness in M x N independent variables for 
a grid dimension of MxN. An example of the almost per- 
fect registration of data onto the grid is shown in Fig. 1. 
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Figure 1. Pointing directions observed towards Ceres at 
105 ixm. The spacecraft raster dimension was 3 x 3 in 
the spacecraft coordinates Y x Z. Sky sampling in Y is 
achieved through a combination of chopper scans in Y and 
repointings of the spacecraft. In Z the sky sampling is de- 
termined alone by the spacecraft repointing interval. Each 
pointing direction for 11833 individual data samples for 
the central pixel of the CI 00 detector array is plotted as a 
cross. The rectangular grid is the "P32 natural grid" (see 
text) for this observation, on which the sky brightness dis- 
tribution is to be solved. The grid sampling is 14.99 x 45.98 
arcsec. Only points with the on target flag set (i.e. not 
including slews) are plotted. The pointings typically lie 
within 1 arcsec of the centre of each pixel of the P32 nat- 
ural grid. 

All the example maps shown in this paper will be sam- 
pled without a gridding function on the P32 natural grid. 
This means that the data in each of the map pixels is 
independent. This is unlike the situation for most other 
maps calculated using PI A, where a gridding function is 
employed and so raw data can contribute to more than 
one map pixel. 
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3. TRANSIENT RESPONSE OF THE PHT-C DETECTORSup of charge on the contacts between semiconductor and 



Gallium doped germanium pliotoconductor detectors ex- 
hibit a number of effects under low background condi- 
tions which create severe calibration uncertainties. The 
ISOPHOT C200 (Ge:Ga stressed) and especially the ClOO 
(Ge:Ga) detectors have a complex non-linear response as 
a function of illumination history on timescales of ~ 0.1 — 
100 sec, which depends on the absolute illumination as well 
as the changes in illumination. The following behaviour is 
seen in response to illumination steps: 

— a rapid jump in signal, which however undershoots the 
asymptotic equilibrium response (often severely for the 
ClOO array) 

— a hook response (overshooting or undershooting on in- 
termediate timescales of ^ 1 — 10 s) 

— a slow convergence to the asymptotic equilibrium re- 
sponse of the detector on timescales of ^ 10 — 100 s 

An example showing these effects for a staring observa- 
tion is given in Fig. 2. Also seen are points discarded by 
the deglitching algorithm in PIA, marked by stars. Many 
glitches are accompanied by a longer lived "tail". It is 
these tails that determine the ultimate sensitivity of the 
ISOPHOT ClOO array. 
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Figure 2. Example of the transient response of a Ge:Ga 
detector following an upwards step to a constant level of 
illumination at 100 fim.. The hook response is shown in the 
inset. 



3.1. Semi-empirical detector model 

Unlike the ISOCAM and ISOPHOT-S Si:Ge detectors, the 

transient response behaviour of the Ge:Ga detectors of 
ISOPHOT-C could not be described by the Fouks Schu- 
bert model (Schubert et al. 1995; Acosta-Pulido 1998). 
However, the effects observed in ISOPHOT-C can to some 
extent be qualitatively understood in terms of the build 



readout circuit, as modelled by Sclar (1984). Sclar's model 
can be formulated as the superposition of three exponen- 
tials to describe the response of the detector to a step 
in illumination. The three timescales involved are a mini- 
mum requirement to reproduce the hook response and the 
long term approach to the asymptotic response. However, 
no theory exists to predict values for the constants in this 
representation, and thus they must be found empirically. 

In practice it is a formidable and unsolved problem 
to extract the constants from in flight data for as many 
as three temporal components to the transient response 
behaviour. In this paper we describe an approximate so- 
lution involving a superposition of two exponentials, each 
with illumination-dependent constants. This was found to 
adequately represent the detector response on timescales 
greater than a few seconds, though it only approximately 
models the hook response. Our model can be viewed as a 
generalisation of the single exponential model developed 
for correction of longer stares by Acosta-Pulido, Gabriel 
& Castaneda (2000). 

In the model considered here, the basic detector re- 
sponse S to an instantaneous step in illumination from 
5oop to Soo is given by the sum of a "slow" signal response 

and a "fast" signal response ^2 as: 

^ - 5i + ^2 (1) 
Si = (l-/3i)^oo(l-cxp[-t/Ti]) + ^01 exp[-t/ri] (2) 

^2 = 132 Sac (1 - exp[-i/T2]) + So2 exp[-i/r2] (3) 

Soi and S02 are signals corresponding to the slow and fast 
response components immediately following the illumina- 
tion step. Soi and 502 are related to the corresponding sig- 
nals immediately prior to the illumination step. Sip and 
S'2p, which contain the information about the illumination 
history. 5*01 and 5*02 are specified by the jump conditions: 



Soi = Pi {Soo — Soop) + Sip (4) 
'S'02 — S2p (5) 

The model is thus specified by four primary parameters: 
the slow and fast timescales ri and T2, the "jump" factor 
(3i and a constant of proportionality for the fast response 
(32- In order to obtain satisfactory model fits to complex 
illumination histories, each of these primary parameters 
had to be specified as monotonic functions of illumination, 
each with three subsiduary parameters: 

f3i = Pio + fill * S!^' (6) 

Tl = no + Til * (7) 
02 = 020 + 021 * S^^ (8) 



T2 



T20 + T21 * ST, 



(9) 



Thus 12 subsiduary parameters in all are needed to char- 
acterise the transient response behaviour of the detectors. 
In addition, the initial starting state of the detector, given 
by the values of Sip and 52p immediately prior to the start 
of the mapping observation, must be specified. 



4 



3.2. Determination of model parameters 

The paramctcrisation of the illumination dependence of 
the primary parameters ri , /3i , T2 and /32 was made in en- 
gineering units (V/s), rather than in astrophysical units 
such as MJy/sr. Thus, it was assumed that there is no 
dependence of the detector response on the wavelength of 
the IR photons reaching the detector. Since the individual 
pixels of the CI 00 and the C200 arrays behave as indepen- 
dent detectors, each of the 12 parameters of the detector 
model had to be separately determined for each detector 
pixel. 

The "slow" primary parameters ti and Pi were found 
from fits to the detector response to fine calibration source 
(FCS) exposures. An example is given in Fig. 3. The model 
fit to the signal response on time scales ranging from a few 
seconds up to the duration of the FCS exposure depends 
mainly on the values of ti and A database of val- 
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Figures. Example of a model fit (curve) to the response 
of the central pixel of the CI 00 array to a constant FCS 
illumination starting at t = s. Prior to the FCS illumi- 
nation the detector was viewing blank sky. The horizontal 
line indicates the fitted value of the illumination. No data 
was recorded in the first two seconds following FCS switch- 
on. 
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Figure 4- Values of the parameter f3i for the central pixel 
of the CI 00 array found from fits to the detector response 
to FCS exposures, plotted versus FCS illumination. The 
solid line is a fit of the function f3i ~ /3io + PuS^'^ with 
/3io = 2.12, Pii = -1.82, and (3i2 ^ 0.022. These are the 
parameters adopted for th model giving the default illumi- 
nation dependence of /3i for this pixel. 
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Figure 5. As Fig. 4, but for the parameter ti . The solid line 
is a fit of the function ri = rio -|- tuS'^^ with rio = 6.92s, 
ni = 4.28, and T12 = —1.22. These are the parameters 
adopted for the model giving the default illumination de- 
pendence of Ti for this pixel. 



ues for Ti and /3i was built up for a wide variety of FCS 
illuminations, allowing the dependence of ri and /3i on il- 
lumination to be determined according to Eqns. 6 & 7. An 
example showing the illumination dependence of /?i for 
the central pixel of the ClOO array is given in Fig. 4. The 
solid line is the fitted illumination dependence of /3i used 
to find the values of the parameters /3io, /3ii and /?i2 used 
in the model. A corresponding plot for the illumination 
dependence of n is given in Fig. 5. At shorter timescales 
the response is primarily determined by the "fast" pri- 
mary parameters T2 and ^2- However, the FCS exposures 
could not be used to find T2 and P2 since no data was 
recorded in the first two seconds following the detector 
switch-on. This is also the reason why the hook response 



is not seen in Fig. 3. Furthermore, there is very proba- 
bly a non-negligible timescale needed by the FCS to reach 
its final operating temperature. Therefore, the "fast" pa- 
rameters (with timescales ~0.1< r <~5s) were found 
using a self calibration technique (described in Sect. 4.2), 
operating on co-added repeated chopper sweeps crossing 
standard point source calibrators. Observations of stan- 
dard point source calibrators of a range of flux densities 
allowed the dependence of ti and /?i on illumination to 
be determined according to the parameterisation of Eqns. 
8 & 9. Examples showing the illumination dependence of 
/?2 and T2 for the central pixel of the ClOO array are given 
in Figs. 6 and 7, respectively. 



5 



1.0 7 
0.8 - 



0.6 - 




0.01 0.10 1.00 10.00 100.00 

ILLUMINATION (V/s) 

Figure 6. Values of the parameter for the central pixel 
of the ClOO array found from self calibration optimisation 
of the model on observations of standard point source cal- 
ibrators, plotted versus illumination. The solid line is a fit 
of the function (32 = ^20 + fhiS^l^^^ with P20 = -0.543, 
P21 — 0.723, and (322 = —0.0 1 03, where SiHum is the known 
illumination of the detector pixel due to the point source 
calibrator. These are the parameters adopted for the model 
giving the default illumination dependence of (32 for this 
pixel. 
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Figure 7. As Fig. 4; but for the parameter T2. The solid 
line is a fit of the function T2 = T20 + T2iSi,^^^ with 
T20 = 14.89 s, T21 = -14.24, and T22 = 0.01025. These 
are the parameters adopted for the model giving the de- 
fault illumination dependence of T2 for this pixel. 

4. THE P32TOOLS ALGORITHM 

In a standard reduction of ISOPHOT data using the in- 
teractive analysis procedures of PIA, the analysis is made 
in several irreversible steps, starting with input "edited 
raw data" at the fuh time resolution, and finishing with 
the final calibrated map. To correct for responsivity drift 
effects, however, it is necessary to iterate between a sky 
map and the input data at full time resolution. This dif- 
ferent concept required a completely new data reduction 
package. In its development phase, and for the determina- 



tion of the detector model parameters, this data reduction 
package was run as a set of IDL scripts. It was later in- 
terfaced to PIA as P32TOOLS using a GUI interface as 
described elsewhere in this volume. 

There are three basic elements to the P32TOOLS con- 
cept. The first step is signal conditioning, which aims to 
provide a stream of signal values, each with an attached 
time, from which artifacts such as glitches, dark current, 
and non-linearity effects have been removed, but which 
retains the full imprint of the transient response of the 
detector to the illumination history at the sharpest avail- 
able time resolution. We will refer to such a stream of 
signals as the "signal timeline" . 

The second step is the transient correction itself, which 
is an iterative process to determine the most likely sky 
brightness distribution giving rise to the observed signals. 
This is a non-linear optimisation problem with the val- 
ues of sky brightness as variables. For large maps several 
hundred independent variables are involved. 

Lastly, there is the calculation of calibrated transient 
corrected maps. This entails the conversion from V/s to 
MJy/sr, flat fielding, and the coaddition of data from 
different detectors pixels. Within P32TOOLS, this third 
stage is performed using standard PIA software, and is 
therefore not described further here. 

4.1. Signal Conditioning 
The processing steps for signal conditioning are as follows: 

1. Edited Raw Data is imported, after having been cor- 
rected in PIA for the integration ramp non-linearity. 

2. The specification of chopper phase in the data is checked 
and corrected where necessary. 

3. The pointing is established from the instantaneous satel- 
lite pointing information, unlike PIA which uses the 
commanded pointing. This also establishes the point- 
ing directions on the slews between the spacecraft fine 
pointings. Internally defined "on target flags" are set 
for each spacecraft fine pointing direction, and the cen- 
tral direction and sampling intervals directions for the 
P32 "natural grid" (see Sect. 2) are established. 

4. The integration ramps are differentiated and the data 
is corrected for dark current and dependence of pho- 
tometry on reset interval, using standard calibration 
parameters taken from PIA. 

5. Systematic variations of the signal according to po- 
sition in the integration ramps are removed from the 
data. This process is called "signal relativisation" within 
P32TOOLS. In contrast to the default pipeline or PIA 
analysis, this allows the use of all non-destructive reads, 
as well as providing a more complete signal timeline 
for use in the correction of the transient behaviour of 
the detector. The "signal relativisation" is applied such 
that signal variations due to the transient response be- 
haviour are preserved. 
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Figures. Detail from a signal timeline from an observation of the standard calibrator Ceres, observed in the C105 
filter. The green line shows the fit to the observed signal (shown in red) from pixel 1 of the ClOO array, as found 
by the P32TOOLS algorithm for correction of the transient response behaviour of this pixel. 5 chopper sweeps, each 
comprising 13 pointing directions ("chopper plateaus") are shown, divided by vertical dotted lines separated in time by 
a duration of 6.1s for each sweep. The upper line in the figure shows the fine pointing flag; whereas the first chopper 
sweep shown was made while the .spacecraft was slewing between fine pointings (flag value 0), the subsequent chopper 
sweeps were made on a single fine pointing (flag value 1). The remaining horizontal lines denote masks at the full 
time resolution of the data denoting glitches, readout status and whether or not a solution for the illumination could 
be found by the algorithm. The overall solution for the sky illumination for this detector pixel is given by the histogram 
in black. As described in the text, this overall solution is calculated from a combination of individual solutions for 
illumination on each successive chopper plateau, which are shown here as purple bars. For the very brightest sources, 
as shown in this example, the derived illuminations can be up to a factor of 6 brighter than the raw data for pixels 
in the ClOO array. In this example the hook response is well seen in each chopper sweep at the 9th and 10th chopper 
plateaus. 



4.2. Transient Correction 



The calculation to correct for the transient response be- 
haviour is made separately for each detector pixel. The 
kernel of the procedure is to solve for the illumination 
corresponding to the signal on each chopper plateau. This 
is done in a binary chop scheme for set values of the 12 



Random noise is determined by examining the statis- 
tics of the signals on each chopper plateau. This is nec- 
essary due to the high readout rates of the detectors in 
many P32 observations. Otherwise, the determination 
of noise from individual readout ramps, as done in a 
standard PIA analysis, would be insufhcicntly precise. 
Spikes are detected and removed from the data in a 

first-stage deglitching procedure. This operates at the , ^ , mi i r ^ i ^ i i 

„ , . , ^. i-ii 1 i • 1 ii 1 i i- detector parameters, ihe values tor bir, and needed 

full time resolution 01 the data given by the non-dcstructive , r .,, . . , , , , 

to solve for illumination on each chopper plateau are cal- 
culated from the response of the model to the preceding 
illumination history. After solving for the illumination on 
the current chopper plateau, the found value is combined 
with previous solutions encountered for the same sky po- 
sition from previous chopper plateaus viewing the same 
sky direction. The response of the detector pixel to the il- 
lumination history is recalculated and a solution is found 
for the illumination viewed on the next chopper plateau. 



read interval. 

A second stage deglitching procedure is performed, op- 
erating at the time resolution of the chopper plateau 
(which comprise of at least 16 non-destructive read in- 
tervals) . This procedure (described by Peschke & Tuffs 
in this volume) removes any long lived glitch "tails" 
following the spikes detected in the first-stage deglitch. 
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Figure 9. Detail from a timeline of signal versus time from another observation of Ceres, this time observed in the 
C160 filter using the C200 detector array. The key to the plot is as given in the caption to Fig. 8. This example 
shows the transition from a spacecraft fine pointing for which the chopper is sweeping through beam sidelobes to a fine 
pointing where the chopper samples the beam kernel. In this case the red line corresponding to the data is not seen 
during the intervening spacecraft slew, as it is exactly overplotted by the green line showing the fit to the data (see 
text). The chopper sweeps on the second fine pointing encompassing the beam kernel show a typical enhancement of 
the source-background contrast induced by the algorithm for correction of the transient response. 



This procedure is continued until all the signal timeline 
has been processed. Then a second pass through the data 
can be made, and so on, until the overall solution for the 
sky brightness distribution no longer changes. The good- 
ness of the solution is quantified through a value for x^, 
calculated from a comparison between the model fit and 
the signal timeline. The random uncertainties assigned to 
each signal sample in the timeline can optionally be re- 
calculated from the self consistency of solutions for the 
illumination found in a given sky direction. Examples of 
model fits to the signal timeline are shown in Fig. 8 for 
the ClOO detector and in Fig. 9 for the C200 detector. 

The algorithm can also be used for datasets for which 
the readout rate was too rapid to allow all the data to be 
transmitted to the ground. In such cases the illumination 
history can still be determined at all times on fine point- 
ings, even where there are gaps in the signal timeline. On 
slews, the illumination history during the gaps are either 
interpolated in time, or taken from the nearest previously 
solved direction on the P32 natural grid. 

For most observations using the C200 detector, and ob- 
servations of moderate intensity sources made using the 
ClOO detector, good fits to the entire unbroken signal 
timeline can be obtained for each detector pixel. For obser- 



vations of bright sources using the ClOO detector, however, 
it is often necessary to break up the solutions into smaller 
segments, according to the distribution of the sources in 
the field mapped. This is best achieved by solving sepa- 
rately for the portions of the signal timeline passing through 
the bright sources. The most conservative processing strat- 
egy, appropriate for observations using the ClOO detector 
of complex brightness distributions with high source to 
background contrast, is to break up the calculation into 
time intervals corresponding to individual spacecraft fine 
pointings. This means that no processing of the slews is 
required. 

Several processing options are available, to optimise 
solutions for the particular characteristics of each dataset: 

Determination of detector starting state 

The default operation of the program is to assume 
that the detector is in equilibrium at the start of the 
observation, viewing the sky direction corresponding to 
the first detector plateau. This corresponds to the case 
Sip = Soi = Soo and S'lp = 0. This will very rarely be 
the case, however, since an FCS measurement which may 
not be perfectly matched to the sky brightness will have 
been made just seconds prior to the start of the spacecraft 
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raster. Therefore, the program can optionally find a solu- 
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Figure 10. A further detail from the timeline of signal ver- 
sus time for the observation of Ceres in the CI 60 filter. 
This time the fit to the data for the very first fine pointing 
is shown, to illustrate the fit to the .switch-on transient 
behaviour of the detector (see Sect. ^.2) For a constant 
repeated illumination pattern in each chopper sweep, both 
the data (red) and the fit (green) gradually move lower 
through the first half minute of the observation. 

tion for the detector starting state, as characterised by the 
values of Sip and 5'2p immediately prior to the start of the 
first chopper plateau. This is done by solving for Sip and 
iS'2p for a fixed 6*00 on the first plateau, where Soo is deter- 
mined from the last chopper plateau in the first spacecraft 
pointing to view the same sky direction. The timeline for 
the first spacecraft pointing is then repeatedly processed 
until a stable solution for S'lp and S'2p prior to the first 
plateau is found. This option was found particularly use- 
ful for data taken with the C200 detector. An example of 
a model fit to a timeline where a solution for the detector 
starting states was found in this way is shown in Fig. 10. 

"Self calibration" of detector parameters 

In order to determine the primary detector parameters 
(32 and T2 in the first instance, "a self calibration" process- 
ing technique was used by which any of the 12 detector 
parameters (Sect. 3.1) could be determined from observa- 
tions of bright sources. Sources of any arbitrary brightness 
distribution and position can be used. However, to deter- 
mine the detector parameters for the "fast" response of 
the detector, point source calibrators of known brightness 
were used, and the solutions for T2 and (32 were addition- 
ally constrained by the requirement that the solutions for 
illumination were consistent with the known flux densities 
of the calibrators. 

The self calibration procedure works by simply repeat- 
ing the optimisation process for different subsets of the 12 



detector parameters, specified by the user. Each combina- 
tion of detector parameters produces an individual solu- 
tion for the sky brightness distribution, as well as a value 
for . A search is made in the parameter space of the 
detector parameters until a minimum is found in x^. This 
is a lengthy process, however, making it advisable to per- 
form self calibrations on limited portions of the timeline. 
This is often chosen to correspond to a single spacecraft 
pointing where the chopper sweeps for the detector pixel 
being investigated pass through bright structure. The pro- 
cess can be further accelerated in cases that the repeated 
chopper sweeps through a source give the same repeated 
signal pattern. Then, the optimisation can be performed 
on an average of the coadded chopper sweeps performed 
on a single spacecraft pointing. This latter technique is re- 
ferred to in P32TOOLS as a "composite self calibration" . 
Some examples of the useage of the self calibration proce- 
dure are given in this volume by Schulz et al. (2002). 

Processing of slews 

The calculation of the illumination history for data 
taken during the slews is less straightforward than for 
data taken during the spacecraft fine pointings. This is be- 
cause the instantaneous pointing directions during slews 
are effectively at random positions with respect to the 
P32 natural grid. Therefore, the use of solutions for the 
sky brightness distribution on the P32 natural grid found 
from analysis of previous fine pointings will in general lead 
to inaccuracies in the determination of the illumination 
history during the slews. There are two main approaches 
which can be adopted to counter this problem. 

The first is to solve for the illumination not over a 
chopper plateau, but for individual readouts in the signal 
timeline during the slew. This automatically provides a 
solution in which the model fit exactly matches the data, 
as there is only one fitted data point for each determined 
illumination. An example is given in the fit to the sig- 
nal timeline for the C200 detector shown in Fig 9. This 
approach is useful for observations of bright structured 
sources in which the variations in illumination with posi- 
tion on the slew on timescales of a chopper plateau exceed 
the signal to noise on individual data samples. It is the de- 
fault processing option for the C200 detector. 

In cases of low signal to noise, and for all instances 
where not all data are transmitted to the ground, solving 
for the average illumination on each chopper plateau is 
a preferable technique. Otherwise, the values of S'lp and 
S'2p calculated just prior to the first chopper plateau on 
the next fine pointing can be wildly inaccurate, which can 
corrupt the solutions for the remainder of the signal time- 
line of the observation. Processing on a plateau by plateau 
basis is the default processing option to calculate the il- 
lumination history during slews for observations with the 
CI 00 detector. 
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5. PHOTOMETRIC PERFORMANCE 

In general the corrections in integrated flux densities made 
by the algorithm depend on the source brightness, struc- 
ture, the source/background ratio, and the dwell time 
on each chopper plateau. The largest corrections are for 
bright point sources on faint backgrounds. 

Here we give as an example results achieved for the 
faint standard star HR 1654. This source was not used 
in the determination of the detector model parameters, 
so it constitutes a test of the photometric performance of 
ISOPHOT in its P32 observing mode. The derived inte- 
grated flux densities, with and without correction for the 
transient response behaviour of the CI 00 detector, were 
compared in Fig. 11 with predicted flux densities from a 
stellar model. The corrected photometry is in reasonable 
agreement with the theoretical predictions. As expected, 
there is a trend for observations with larger detector il- 
luminations to have larger corrections in integrated pho- 
tometry. 
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Figure 11. The measured integrated flux densities of the 
faint standard star HR 1654 plotted against the pre- 
dicted flux densities from a stellar model. The observa- 
tions were done in various filters using the CI 00 detector. 
The photometry derived with and without processing with 
P32TOOLS is shown with stars and crosses, respectively. 



A good linear correlation is also seen between inte- 
grated flux densities of Virgo cluster galaxies (Tuffs et 
al. 2002), derived from P32 ISOPHOT observations pro- 
cessed using the P32TOOLS algorithm, and flux densi- 
ties from the IRAS survey (Fig. 12). The interpretation of 
these measurements constitutes the first science applica- 
tion (Popescu et al. 2002) of P32TOOLS algorithm. The 
ISOPHOT observations of Virgo cluster galaxies were fur- 



thermore used to derive the ratios of fluxes measured by 
ISO to those measured by IRAS. The ISO/IRAS ratios 
were found to be 0.95 and 0.82 at 60 and 100 /im, respec- 
tively, after scaling the ISOPHOT measurements onto the 
COBE-DIRBE flux scale (Tuffs et al. 2002). 
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Figure 12. Integrated and colour- corrected flux densities of 
Virgo cluster galaxies m,easured in the ISO CI 00 filter ver- 
sus the corresponding flux densities measured by IRAS in 
its 100 iim band (taken from Fig. 7 of Tuffs et al. 2002). 
The dotted line represents the relation ISO/IRAS=0.82. 
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Figure 13. Brightness profiles along the Y spacecraft di- 
rection through the standard star HR 1654 o,t lOO^m. 
The solid line represents the brightness profile obtained 
after processing with P32TOOLS, while the dotted line 
shows the profile derived from identically processed data, 
but without correction for the transient response of the de- 
tector. 

Fig. 13 shows an example of a brightness profile through 

HR 1654 at 100 /im, for data processed with and without 
the responsivity correction. Some 95% of the flux density 
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has been recovered by P32TOOLS. Without the correc- 
tion, some 30% of the integrated emission is missing and 
the signal only reaches 50% of the peak illumination. The 
local minimum near 60 arcscc in the Y offset is a typical 
hook response artifact, where the algorithm has overshot 
the true solution. This happens for rapid chopper sweeps 
passing through the beam kernel. This is a fundamental 
limitation of the detector model, which, as described in 
Sect. 3.1, does not correctly reproduce the hook response 
on timescales of up to a few seconds. This problem is par- 
ticularly apparent for downwards illumination steps. The 
only effective antidote is to mask the solution immediately 
following a transition through a bright source peak. An- 
other effect of the inability to model the hook response is 
that the beam profile becomes somewhat distorted. This 
also has the consequence, that for observations of bright 
sources, the measured FWHM can become narrower than 
for the true point spread function, as predicted from the 
telescope optics and pixel footprint. An example of an ex- 
tremely bright point source showing this effect is given in 
Figs. 14 & 15, depicting maps of Ceres in the C105 filter, 
respectively made without and with the correction for the 
transient response of the ClOO detector. 
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Figure 15. Contour m,ap of Ceres in the CI 05 filter with- 
out correction for the transient response behaviour of the 
detector. 50 linear contours have been plotted between 
11.5 and 334 MJy/sr. Measured integrated flux density 
after background subtraction is 28.0 Jy± 0.07 Jy (random) 
±4% (systematic). 



of resolved sources. Also visible in the corrected map is a 
trace of a beam sidelobe. 
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Figure 16. The interacting galaxy pair kpg 347 observed in 
the C200 filter after processing with P32TOOLS. 



Figure 14- Contour map of Ceres in the C105 filter af- 
ter responsivity drift correction. 50 linear contours have 
been plotted between 33 and 2519 MJy/sr. Measured 
integrated flux density after background subtraction is 
98.4: Jy ±0.2 Jy (random) ±6% (systematic). The actual 
flux density of this standard calibrator (from a stellar 
model) is 109 Jy. 

Despite the limitations due to the lack of a proper mod- 
elling of the hook response, the algorithm can effectively 
correct for artifacts associated with the transient response 
on timescales from a few seconds to a few minutes. This is 
illustrated in Figs. 16 & 17 by the maps of the interacting 
galaxy pair KPG 347 in the C200 filter (again, after and 
before correction for the transient response of the detec- 
tor, respectively). The uncorrected map shows a spurious 
elongation in the direction of the spacecraft Y coordinate, 
which is almost completely absent in the corrected map. If 
uncorrected, such artifacts could lead to false conclusions 
about the brightness of FIR emission in the outer regions 



The image of M 101 in Fig. 18, made using the ClOO 
detector with the ClOO filter, is given as a state of the art 
example of what can be achieved with a careful interactive 
processing of P32 data using P32TOOLS. In addition to 
the transient response corrections and a masking of resid- 
ual hook response artifacts, a time dependent flat fleld has 
been applied. The spiral structure of the galaxy, with em- 
bedded HII region complexes and a component of diffuse 
interarm emission can clearly be seen. 
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Figure n. The interacting galaxy pair kpg 347 with identi- 
cal processing, except that the correction for the transient 
response behaviour of the detector has been omitted 
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